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Abstract 



In this paper we discuss bimeron pseudo spin textures for double layer 
quantum hall systems with filling factor u = 1. Bimerons are excitations 
corresponding to bound pairs of merons and anti-merons. Bimeron solutions 
have already been studied at great length by other groups by minimising the 
microsopic Hamiltonian between microscopic trial wavefunctions. Here we 
calculate them by numerically solving coupled nonlinear partial differential 
equations arising from extremisation of the effective action for pseudospin 
textures. We also calculate the different contributions to the energy of our 
bimerons , coming from pseudospin stiffness, capacitance and coulomb in- 
teractions between the merons. Apart from augmenting earlier results , this 
allows us to check how good an approximation it is to think of the bimeron 
as a pair of rigid objects (merons) with logarithmically growing energy, and 
with electric charge ^. Our differential equation approach also allows us to 
study the dependence of the spin texture as a function of the distance be- 
tween merons, and the inter layer distance. Lastly, the technical problem of 
solving coupled nonlinear partial differential equations, subject to the special 
boundary conditions of bimerons is interesting in its own right. 
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I. INTRODUCTION 



The fascinating discoveries of the quantum Hall effect (QHE) originally found in single 
two dimensional electron layers, have also been extended to double layer systems , thanks to 
the development of techniques for growing GaAs heterostructures containing two separated 
layers of two-dimensional electron gas (see for example references Apart from finding 
plateaus in Hall conductivity at total filling fractions v corresponding to the "direct sum" 
of the familiar integral and odd-denominator fractional QHE in the individual single layers, 
experiments also show the occurence of new plateaus which are intrinsic to double-layer 
systems and rely on interlayer quantum coherence and correlations. On the theoretical 
front, a large body of work has already been done on double-layer systems. An extensive list 
of references to this literature has been given in the lucid review of this subject by Girvin 
and MacDonald and in the paper by Moon etal P]. 

Generally one analyses double layer systems by attributing to the electrons, in addition 
to their spatial coordinates on the plane, a two-component "pseudospin" whose up and 
down components refer to the amplitude for the electron to be in the first and second layers 
respectively. The real physical spin of the electrons is assumed , as a staring approximation, 
to be fully polarised by the strong magnetic field and hence frozen as a degree of freedom. 
However, even when real physical spin is suppressed, the use of a pseudospin to represent the 
layer degree of freedom maps the double layer spinless problem into a monolayer problem 
with spin @]. Such a mapping allows one to borrow for double layer systems, the rich body 
of insights and results available from single layer systems with real spin. 

Thus one may expect a fully symmetric (polarised) pseudospin state to be energetically 
preferred because of a combination of Coulomb repulsion and the Pauli principle which 
forces an associated antisymmetric spatial wavefunction, just as in itenerant ferromagnetism. 
Further, the relevance of Skyrmions to systems with real spin, predicted by theoretical 
considerations and supported by experimental evidence 0, has in turn prompted 
studies of similar topological excitations in spinless double layer systems, but now involving 
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pseudospin (See Girvin and MacDonald p|, Moon etal and references given therein). 

Because of interplane-intraplane anisotropy in Coulomb repulsion between electrons lo- 
cated in the two layers, as well the capacitance energy of maintaining unequal charge den- 
sity in the two layers, the effective Action governing pseudospin enjoys only U(l) symmetry 
of rotations about the z-axis (the direction perpendicular to the x-y plane of the layers). 
Finiteness of the capacitance energy between the two layers requires that asymptotically 
the pseudospin must lie on the easy (x-y) plane. The basic topological excitations in that 
case are the so-called merons which are vortices in pseudospin with a winding number of 
one-half (with respect to the second Homotopy group 112. These are similar to vortices in 
the X-Y model, but non singular at the origin since the pseudospin is not restricted to lie 
on the x-y plane. But like the former they do have an energy that grows logarithmically 
with size. One can also have meron anti-meron bound pairs whose energy is finite. Such 
a pair is topologically equivalent to Skyrmions and carries unit winding number. (For an 
introduction to such topological excitations, their winding numbers, etc. see reference [§.) 

The possibilty of topological excitations like merons and bimerons in double layer systems 
has generated much interest, in part because of the excitement surrounding the Skyrmion 
excitations in systems with real spin, and in part because of the additional possibility here 
of Kosterlitz-Thoulles type phase transitions caused by the break-up of bound bimerons 
into separated meron pairs [0,0- Bimeron solutions have already been extensively studied 
in a body of papers by Girvin, MacDonald and co-workers ]ri| and . These calcu- 
lations are based on optimising microscopic wavefunctions with respect to the microscopic 
interaction Hamiltonian. 

We will also calculate bimeron solutions and their energies here, but by using an alter- 
nate method. An Effective Action for slowly varying pseudospin textures has already been 
obtained by Moon et al [Q]. If one extremises that Action one will get differential equations 
which the unit-vector valued field of psuedo spin configurations m(r) should obey in the 
classical limit. 

In this paper we solve these coupled non-linear differential equations, through a com- 
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bination of analytically motivated ansatze followed by numerical calculations. We obtain 
bimerons as approximate time-independent solutions with appropriate topologically non- 
trivial boundary conditions, for a range of separations between the meron and its partner 
the anti-meron and also for a set of different inter-layer distances. The dependence of the 
bimeron texture on these variables is discussed. They turn out to be reasonably similar to 
what one would expect on general grounds. We also obtain the energy of this bimeron as a 
function of the separation between the meron centers. We include in this energy contribu- 
tions coming from the pseudospin stiffness, its anisotropy, the capacitance energy and the 
Coulomb energy. By minimising this energy with respect to the meron separation, we are 
also able to give an independent value for the optimal meron separation in a bimeron. We 
compare these results with earlier work, including our own. 

Apart from this, our work also enables us to independently check the validity of a physical 
picture often used |]ll| in estimating bimeron energies, namely, that they can be viewed as 



a pair of rigid objects carrying electric charge of ^ and a logarithmically growing energy. 
A work somewhat similar in spirit to ours, but in the context of Skyrmions of real spin 
systems was done by Abolfath et.al. who compared results obtained from solving a non- 
linear differential equation with those obtained from microscopic calculations For yet 



another way of approaching meron solutions, starting from a Chern- Simons field theory see 



the work of Ichinose and Sekiguchi [jT3|. 



In an earlier paper |jT^ we had done a similar study of single meron solutions. But the 



present work is much more complicated at the computational level. Single meron solutions 
are circularly symmetric, with the spin component on the plane pointing along the coordinate 
direction. Thus the only unknown, namely, the spin-z component obeys an ordinary (though 
non-linear) differential equation in the radius variable r. Further, the boundary conditions 
relevant to a single meron can be imposed, conveniently, at the end points r = and r = oo. 
By contrast the boundary conditions characterising a bimeron are rriz = ±1 at two finite 
points on the plane where the two merons have their centers. The spin direction is also not 
simply related to the coordinate direction, so that there are two independent fields, say, 



and tan~^y^j, (since the m is constrained to be a unit vector) which obey coupled partial 
differential equations on the plane. We found it quite challenging to analyse these coupled 
equations analytically as far as posible, and use that information to employ an appropriate 
ansatz and coordinate system to numerically solve the equtions (using a desk-top computer). 

Finally, we should reiterate that our work here clearly relies heavily on the advances 
already made by the Indiana group 0, [|T^ and is to be viewed as something which 
will hopefully augment their findings. 



II. THE SPIN TEXTURE EQUATIONS 



The differential equations obeyed by spin textures is obtained by extremising an effective 
action which has already been derived by Moon et al starting from the basic microscopic 



physics. See also Ezawa [I^. These results were summarised in our earlier paper [|14 



Briefly , the pseudospin texture of a state is described by a classical unit vector m{r) which 
gives the local direction of the pseudospin. Here r is the coordinate on the x-y plane 
carrying the layers, while the magnetic field B is along the z-direction. The fully polarised 
"ferromagnetic" ground state corresponds to m pointing everywhere in the same direction, 
say, along the x-axis. Using this as the reference state, any other state with some arbitrary 
texture m{f) is given by performing a local pseudospin rotation on this uniform ground 
state . The leading low-wavelength terms in the effective Action for time independent 
configurations m(r), as obtained by Moon et al 0] is 



J(m) : 
where 



+ Ci[m] + C2[m] (2.1) 



drdr'V {f — r')q{r)q{r' 



(2.2) 



and 



e'^d^ /■ 2 /■ ,2 ,.m^(r)V2m^(r') 
62 m = — / d r / d r ; 



(2.3) 



The constants pA and pe are pseudospin stiffness parameters whose physical origin is the 
exculsion principle (Hund's rule) mentioned earlier. They are given by 



PA = (7^) / dke Vjf exp{—- 
6It( Jo I 



PE = (7777^) / Vi exp{^) (2.4) 



^327r^'yo " ^2 

where 14"^ = 27re^/(eA;) and V^^ = (exp(—A;(i)27re^)/(eA;) are the Fourier transforms of the 
Coulomb interactions between electrons in the same and different layers respectively. All 
distances (and inverse wave vectors) are in units of the magnetic length /. The jSml term 
represents the so-called capacitance or charging energy needed to maintain unequal amounts 
of charge density in the two layers. Recall that the z-component of pseudospin represents 
the difference between the densities in the two layers. The constant /5 is given by 

/?=(—) / dkk {V%0) - V^{k)) expi^) (2.5) 

where Vj^ = ^{Vjf — Vjf). Finally, q{r) is the topological density associated with pseudospin 
texture, which is also its charge density 0. It is given by 

= -^^vf,m{f).[d^in{f)xdf,m{f)] (2.6) 

OTT 

We will refer to the the non-local term Ci, as the Coulomb term since it has been identified 
as the Coulomb energy associated with topological structures in the pseudospin textures , 
0. The other non local term C2 arises in the gradient expansion but is not amenable to 
simple physical interpretation. 

The field equations are obtained by extremising this Hamiltonian with respect to the 
independent field variables, which can be taken to be and a = tan~^ {^^t^ ■ This a is 
just the azimuthal angle of the projection of m on to the x-y plane. The non-local terms Ci 
ansd C2 in the Action ( p.l| ) will render the field equations into coupled integro-differential 
equations. While in the single meron case we did solve such an integro differential equation 
|T^ , for the more complicated case of bimerons we will be content to solve the equations in 



the absence of the integral terms Ci ansd C2 .The contributions of these terms can however 
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be included in the total energy, but by using solutions of the local equations. In mild 
justification of this strategy, we will find later that the Coulomb energy for instance is 
less than half the energy from the local terms in eq. 

The coupled field equations for and a = tan~^ {j^T^ resulting from eq. ( |2.1| ) in the 
absence of Ci and C2 are 

p^V'm, + pEm, J- ^— + + V'a - 2j3m, = (2.7) 

and 

V.[(l -m^)Va] = (2.8) 



III. BIPOLAR COORDINATES 

To find bimeron solutions we have to numerically solve the coupled partial differen- 
tial equations (PDE) in ( p.7|) and ( p.8|) . The defining boundary condition of a bimeron is 



rriz = ±1 at the points (0, ±a) Our strategy will be to use the known exact solution of 
these equations in the Non-Linear Sigma Model (NLSM) limit, and solve the full equations 
iteratively starting with the NLSM solution. The NLSM limit is realised when the layer 
separation d goes to zero in which case we see from their defining equations above that 
Pa = PEi i-e. the stiffness is isotropic and further that the capacitance coefficient (3 van- 
ishes. Then , with Ci and C2 also neglected, the action in ( p.l| ) is just that of the NLSM, all 
of whose solutions are exactly known [Q. They are conveniently described by the complex 
field w(z) which represents the stereographic projection of the unit sphere of textures m. It 
is defined by 

w{z) ^ 7f f 3.1 

(1 - m-^) 

where z = x+iy. Our texture variables and a are related to w(z) by 

\w\^-l 



|W|2 + 1 

and a = arg (w) (3.2) 



Any analytic function w(z) will be a solution of the NLSM. In particular the function 

w{z) = —— (3.3) 
z + a 

represents the bimeron, with the points (0,-a) and (0,a) representing the centers of the two 
merons, where the solution gives = ±1 respectively. It may be checked that ( p.3| ) satisfies 
the coupled equations ( |2.7| ) and ( p.8| ) in the isotropic limit. 

When the interlayer separation d is not zero, we have to cope with the coupled field 
equations ( |2.7D and ( |2.8| ) with both the anisotropic stiffness and capacitance terms present. 
Some analysis of this system was done long ago by Ogilvie and Guralnik [|T^ who studied the 



NLSM with the mass (capacitance) term included but no anisotropy. ( An ansatz suggested 
in ref ( [0]) does not work ,as we will show below.) Meanwhile Watanabe and Otsu |18 



studied the anisotropic NLSM but without the mass term. Both made considerable progress 
analytically, but neither offered exact or numerical solutions. Here we will try to solve ( p.7|) 
and ( |2.8| ) numerically after including both the capacitance and anisotropic terms . 

To do so , it will be convenient to use a bipolar coordinate system to describe the x-y 
plane, as might be expected when we have to impose boudary conditions at two finite points 
(0,-a) and (0,a). These coordinates, rj and 0, are defined by 

.z — a, 

T] = In 



z + a 



and (f) = arg ( ) (3.4) 

\z + aj 

This coordinate set has many advantages |T9[. The points (0,-a) and (0,a) at which we have 



to impose boundary conditions are now mapped into rj ±oo. The full x-y plane is mapped 
in {ri,(f)) coordinates to an infinite strip with rj = [— oo, +oo] and = [— 7r,7r]. Finally, it is 
clear upon comparing eq(|3.4|) to eq ( |3.3| ) that this set of coordinates is closely related to the 
exact NLSM bimeron solution. Clearly the the exact NLSM solution ( |3.3| ) corresponds to 
the simple expressions 

niz = tanhr) 

and a = (j) (3.5) 
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Away from the NLSM limit, since this is an orthogonal coordinate system with simple 
expressions for the gradient, divergence and Laplacian, the equations (|2.7| ) and ( |2.8| ) become 



pE I - ml 



[dim:, + <9!m^) + ^^(^^"^^ "^f't"^^^ + mz{{dr,a + d^af 

2/9 ^2 



Q^(ry,0) = O (3.6) 



(1 — m^)(9^a + 9|a) — 2mz{drjmzdria + d^rrizd^a). = (3.7) 



where 



(coshr/ — cos0j^ 

is the Jacobian of this coordinate transformation. Now let us analyse these equations as 
different terms are included in stages. 

(a) In the NLSM limit, our exact solution has a = (p. Then ( p.7| ) forces to be a 
function of rj alone , tjIz = rriz [rj). Upon inserting this into the other equation ( |3.6|) it 
becomes an ordinary non-linear differential equation . This is the advantage of this choice 
of coordinates. The solution can be verified to be = tanh{ri). 

(b) Next let us include anisotropy (p^ 7^ Pe) , while still keeping the capacitance term 
zero {(3 = 0). Once again we can set a = 0, and consequently = mz{r]), which will obey 
again an ordinary differential equation given by 



, Pa- Pe . ^ 1 



Pe ' I - ml 



ipa N , mz{dr{mzf , „ 

[d + — '—— + mz = (3.9) 

' (1 — m^^Y 

This has no analytic solution, but can be solved relatively easily numerically, being 
just an ordinary differential equation in the variable rj. As boundary conditions we impose 
mz = at = and = 1 at = cxd, (Note that the equation above is symmetric 
under rj —* —t], so that we can choose the solution to be antisymmetric, i.e. m,z{—ri) = 
—mz{ri)). The resulting numerical solutions for different values of layer separation d (on 
which the anisotropy depends), are shown in fig 1. One can see that with increasing the 
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layer separation, and hence increasing anisotropy in the stiffness, the pseudospin component 
rriz reaches its asymptotic value more slowly. 

(c) Finally let us also include the capacitance term and consider the equations ( p.6| and 
|3.7| ) in full. Now the ansatz a = is no longer sustainable, in contrast to what has been 
suggested in ref ( JT^])- The substitution of the ansatz a = in equation ( ^.T] ) would again 
force d^rriz = ,i.e. = mz{r]). But now this is in contradiction with equation ( p.6|) which 
has an explicit dependence through the last (capacitance) term ^ ijli'P)- Therefore 
,once one includes the capacitance term in eg nation ( p. 6|) both a and become functions of 
both 7] and (f). One has unavoidably to solve the coupled non-linear PDE for = rrizii], (p) 
and a = a{rj, 0). 

We do this by employing what we believe is a good ansatz for a which approximately 



satisfies p.7|) . We then solve the other equation (|3.6| ) numerically after inserting that ansatz 
for a. Our ansatz is been motivated by the following arguments. One can see from equa- 
tion( |3.6[ ) that the troublesome dependent term is negligibly small in the large t] region 
{Q ~ sech{r])) and is most dominant in the small rj region. Hence a will still approach as 
1] ^ oo but needs to be modified substantially in the small t] region where however <^ 1. 



When <^ 1 equation ( pj.?] ) can be approximated by 

V^a = (3.10) 

This is just Laplace's equation in two dimension whose solutions are all harmonic functions. 
With this in mind we choose our ansatz for a as follows : 

a = (p — BKexp{—\ri\)sin{(f)) (3-11) 

where 

K = (— )5 a (3.12) 

PE 

This solves Laplace's equation and satisfies all the required boundary conditions and 
asymptotic behaviour, namely 
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±oo 


a = 
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= 


a = 


vr 


when 


(f) = IT 


a = 





when 


K = 



(3.13) 

. Note that the ansatz has a cusp at ?7 = 0. This need not cause concern. Some such 
cusps can be expected on physical grounds and are famihar in sohton physics. The point 
is that each meron feels some force due to the other (Coulomb plus a logarithmic force) 
at arbitrary separation. We would expect them to move because of this force, and cannot 
strictly speaking expect a static bimeron solution to exist at arbitrary separation. But a 
cusp, like the one in the above ansatz, amounts to a delta function in the second derivative 
and can be interpreted as a external force just at rj = which can "hold the two merons 
together" at arbitrary separation. For more discussion of this point see Rajaraman and 
Perring and Skyrme |T5| where this technique was used to get intersoliton forces betwen one 
dimensional solitons. 

The constant B is chosen by minimising the energy. Substituting this ansatz in equation 
i WM we then solved it numerically subject to the boundary condition 

171^=0 at r] = 

= ±1 when i] = ±oo (3-14) 

. It is sufficient to solve the equation in the first quadrant .i.e. (//[O, oo]and(f)[0, tt]) . For the 
rest of the quadrants solutions can be obtained by writing 

a{-V, 4>) = a{v, 0) = -a{v, -0) 

(3.15) 

which is consistent with the invariance of equations ( |3.6| )and( |3.7| ) under the transformation 
1] —7] and — s> —0. 
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IV. NUMERICAL PROCEDURE 



Before proceeding to solve this PDE ( p.6|) we must take note of the fact that the last term 
of the equation ( p^) is singular at the point (r^ = 0, = 0). This point corresponds to spatial 
infinity on the parent x-y plane. As one moves near this point the leading singularity in the 
equation, coming from the term, goes like ,-^2^^2)2 with other subleading singuarities of 
the form —. — ^ . It can be seen that this leading singularity can be offset by requiring that 



m, behave as 



exp 



2k 



. I g(ri,(f)), where this g('r],(f)), is a more smooth function for 

which one solves numerically. This corresponds, in more familiar polar coordinates (r, 6) to 
writing in the form [exp — (^)] g (r, 9). That ruz will suffer such an exponential fall-off 
as r — s> 00 can also be inferred directly from the "mass term " 2(3mz in the original field 
equation ( p.7|) . Similarly one can also verify that the cancellation of the subleading singular 
terms can be achievedd by requiring that g has to behave like \J~ijf + 0^) as 77, 0. 

Given this functional form of near the origin of the 77, plane, the boundary conditions 
( P.14| ), and the ansatz ( |3.11| )for a we solved equation (|3.6|) through an iterative procedure. 
We start with the solution for /t=0 but with full anisotropy, which can be obtained relatively 
easily from the ordinary differential equation ( |3.9| ). We then use this solution as input to 
obtain the solution for k equal to a small number e through the Newton-Raphson method 
|2CI| |. The solution for e is then used as input to obtain the solution for 2e and so on. 
This procedure is repeated until one reaches the desired value of k. The advantage of this 
procedure is that one can make e arbritrarily small to make the Newton-Raphson method 
converge. In this way we obtained solutions for different values of the ansatz parameter B 
for each value of k. 



V. RESULTS AND DISCUSSION 



Our solutions of equations ( p. 6| ) and along with the ansatz ( |3.11| ) give us the value of the 



pseudospin vector m as a function of 77 and 0, or equivalently, the value of the vector-field fh 
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on a lattice of points on the parent x-y plane. We repeated this calculation for a set of values 
of the paramenter B in the ansatz ( |3.11| ). We found that as one varies B starting from 0, 



the energy does not vary much as B goes to 0.1, but then it increases sharply after B=0.1. 
This behaviour is seen to be common to all k and all a. Hence we take B to equal 0.1. and 
solve the PDE for a variety of values of layer separation d, and bimeron separation a (a is 
actually half of the meron-antimeron separation) . Together all these solutions represnt a 
large body of calculated data. But it is neither feasible nor very interesting to try to display 
it all in this paper. Instead we will try to bring out salient features of our solutions through 
examples. 

Recall from ( p.9|) that in the absence of the capacitance term had no - dependence. 
To give some feel for how the varies with in the presence of the capacitance term, we 
plot in fig. 2 the solution m^^r]) of equation ( |3.6| ) for a set of values for 0. This solution 
corresponds to d = 0.7 and a = 3.158. The sequence of curves shown correspond to equal 
to 0, 0.27r, 0.477r, and 0.947r respectively with the outermost one belonging to equal to 0. 
As we have discussed earlier, as rj and tend to zero, the solution should damp exponentially 

as exp( , ). Correspondingly we see in fig. 2 that the low curves rise very slowly as 

7] increases away from zero. We also give for comparison, in the form of the dotted curve, 
the function tanh [r]) which is the solution in the NLSM limit. The comparison shows that 
the restructuring of the pseudospin texture due to the capacitance and term and anisotropy 
is considerable. 

As an alternate representation of our results, we show in fig. 3 the projection of m on the 
x-y plane, for the example of d equal to 0.7 and k, equal to 4.4. (All lengths throughout this 
article are in units of the magnetic length /) . The length of each arrow gives the magnitude 



of its easy-plane projection y + rriy and its direction gives the azimuthal angle of the 
projected vector , namely, a = tan'^ (j^^- "^^^ P^*-*^ clearly shows that our " bimeron" 
solution is indeed a meron-antimeron pair. Note that, as desired, m lies along the x-axis 
asymptotically. This picture closely resembles the general structure obtained by Brey et.al. 
||T0| . The data corresponding to all other values of d and a we studied have a similar behavior. 
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In fig. 4 we plot the topological charge density given in eq.( |2.6| ) as a function of t] and (j) 
in the presence of all the local terms in the field equations, including anisotropic ones. In 
viewing this figure it may be helpful to remember that large |?7| corresponds to the meron 
centers while r/ = 0, = corresponds to spatial infinity. = vr corresponds to the line 
joining the two merons.As topological charge density is symmetric when either of the co- 
ordinate variable changes sign we show the contours only in the first quadrant where both 
1] and are positive. 

Next let us turn to the energetics of these bimeron solutions. In fig. 5. we show how 
the "local" energy i.e. the contribution from the local terms in the energy functional (all 



terms in eq( |2.1| ) except for Ci and C2) varies when one changes the separation 2a between 
the meron and antimeron centres. The appearance of a minimum is quite conspiquous and 
generic to all the layer separations for which the energy is calculated. The example in fig. 5 
corresponds to a layer separation of 0.7. 

In fig. 6 we plot the Coulomb energy Ci evaluated using our solution of the equation (p. 61), 
as a function of the bimeron separation. The continuous curve is the best fit to our calculated 
points . Sometimes in the literature, a phenomenological estimate of bimeron energetics is 
made assuming that it can be viewed as a bound pair of two merons, each symmetrical, 
undistorted by the other and carrying a charge of | . Such a pair would have a Coulomb 
energy of ^ (in units of ^ that we are using). To see how good an approximation this simple 
picture is, we give in the same fig. 6, in the form of a broken line the plot of this function 
^ . We see that the value of the Coulomb energy we get from the actual bimeron solution 
is much larger than what the simple two-charge picture would give. This is presumably 
because each meron is considerably squashed (polarised) by the close proximity of the other. 
In our earlier work on single merons |T^, we had found that at the layer separation {d = 0.7) 



used in fig. 6, the core-radius of individual merons is about 2, which is of the same order 
as the meron-separation in fig 6. In fact we can see that the gap between the two curves 
in fig. 6 is higher for smaller a where the individual merons are squeezed together more. Of 
course our results, while indicative , may not be quantitatively unambiguous. For instance, 
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recall that our solution was obtained using only the local terms in the differential equation 
and the Coulomb energy was calculated by substituting this solution into the integral Ci. 
The non-local Coulomb term's inffuence on the solution has no been included. 

In fig. 7 we plot the variation of three terms in the energy fuctional namely the contri- 
bution from the local terms (capacitance+gradient energy) ,Ci and C2 as a function of the 
bimeron separation. The data presented here corresponds to layer separation d equal to 0.7/ 
but this behaviour is representative of almost all the layer separations ( 0.5, 0.6, 0.7 and 0.8 
for which we have found solutions. The trend of all three contributions is the same for the 
other layer separations also with only slight changes in the slope of the curves. 

Our calculations were done for different bimeron separations a, for each layer separation 
d. In reality, the exact solution should exist only for some optimal bimeron separation a 
for each value of d. One can ask if our calculations would reveal this by minimising the 
total energy at some particular a. To see this, we have shown in Fig. 8 the total energy at 
d=0.7 (i.e. the sum of all three contributions plotted in the fig. 7) as a function of bimeron 
separation a. As we can see from fig. 7, the total energy keeps decreasing with a, all the 
way to about a = 3.2, which is the highest value upto which we could calculate , given 
limitations of our computing facilties. However, the decrease is clearly levelling off and is 
indicative that a minimum may exist at around a=4 or 5. What we have done, in drawing 
fig.8, is to obtain a best-fit-curve of the data points upto a=3.2 and extrapolate that curve 
upto a=4.5. For what it is worth such extrapolation indicates a minimum at about a=4. 
This corresponds to a meron-antimeron separation of of about 8, larger than what Yang 
and MacDonald found by entirely different methods (see their fig. 2) Their value of 

the meron separation for d = 0.7 is about 4.5. We attribute this discrepency to the fact , 
noted already in our discussion of fig. 6, that the Coulomb energy in our explicit calculation 
of the bimeron solution is higher than the undistorted meron pair estimate used in ref( [0). 
The actual larger Coulomb repulsion is, we believe responsible for the larger optimal meron 
separation that we get. 

We saw that the Coulomb interaction energy between the two merons as given by the 
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term Ci in the present calculation differs quite a bit from the simple picture of the bimeron 
as a pair of undistorted merons of charge | each. One can ask if there is a similar discrepency 
in the non-Coulombic energy as well. This is the subject of Table 1. In the picture of a 



bimeron as a pair of merons , |Tl|, |[T^, it will have energy equal to 



Eprev = '2Emc + ^TipE ^'^("j— ) (5-1) 

where Emc and Rmc are respectively the core energy and radius of a single merons, which 
have a logarithmic interaction with each other because of the logarithmic divergence of the 
self energy of single merons. (As stated already we are leaving out their Coulomb interaction 
in the comparison being done in this table.) This E^^-ev has been calculated in our previous 



work It can be compared with the local part of the energy in the present calculation. 
Such a comparison is given in Table 1 for different values of (i, using the optimal value 
of the meron separation a which minimises Eiocai- We see that the comparison is not bad 
considering the completely different ways of estimating this energy in this paper and in 
earlier literature. 

In conclusion, our solution for the bimeron obtained by directly solving the coupled 
partial differential equations that the bimeron texture obeys provides an alternate way of 
obtaining the profiles and energies of these objects. As far as the local part of the energy 
is concerned, the results are in broad agremment with microscopic derivations earlier. But 
the Coulomb energy we obtain is higher by a factor of about 2 from earlier simple estimates 
because in actuality, the two merons in close proximity will not behave like undistorted 
symmetrical merons. 
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FIGURES 

FIG. 1. The solution mz{r]) of equation (|3.9|) . The three continuous curves correspond, as you 
go outwards, to three different values of layer separation d equal to 0.5, 0.6 and 0.7 respectively 
in the unit of magnetic length /. The dotted curve corresponds to the exact solution of NLSM i.e. 
niz = tanh{r]). 

FIG. 2. The solution mz{rj) of equation ( |3.6| ) for a set of values for 0.The curves correspond, 
as you go inwards, to = 0, 0.27r, 0.477r, 0.947r respectively with the outermost one corresponds to 
(j) equal to O.The layer separation d is equal to 0.7Z and bimeron separation a is equal to 3.158Z. 
The dotted curve at the top again corresponds to = tanh{rj). 

FIG. 3. This figure gives the magnitude and direction of x-y projection of m at different points 
on the plane. The layer separation and the bimeron separation are same as in fig. 2. 

FIG. 4. A contour plot of the topological charge density of the bimeron when both the ca- 
pacitace term and the anisotropy term is incorporated. This particular plot corresponds to a layer 
separation dequal to 0.7 and bimeron separation a equal to 3.158 both in the unit of magnetic 
length /. The number against each contour(shown by broken curves) denotes the corresponding 
charge density. 

FIG. 5. This figures gives the plot of the energy (E'/oca/) coming from the local terms in the 
action as a function of bimeron separation a in the unit of magnetic length The unit of energy 

2 

is |j. The points correspond to the actually computed values of the energy while the continuous 
curve is the best fitted curve to it. The form of the best-fit curve \s E = A + B{a — C)^ where A 
and B and C are found out to be .223, .008 and 2.76 respectively. This data corresponds to a layer 
separation d equal to 0.7/ 
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FIG. 6. This figure gives the plot of the coulomb energy as a function of bimeron separation 

2 

a in units of magnetic length /. The unit of energy is jj-. The upper curve is our computed value 
of the coulomb energy integral Ci using the solution of equation ( |3. 61 ) (points). The contious line 
is the best curve to these points. The form of the best fitted curve is E = where A and B are 
found out to be 0.847 and 0.821 respectively. The dotted curve at the bottom corresponds to the 
Coulomb energy that the bimeron would have, if viewed as a bound pair of two point charges of | 
each, separated by a distance 2a . This data corresponds to a layer separation d equal to 0.7 

FIG. 7. this figure gives a relative estimate of the contribution of the three type of terms in 
the action, namely, the local terms, Ci and C2, as a function of bimeron separation a. The units 
are as specified in the earlier figures. This data also corresponds to a layer separation 0.7^ 

FIG. 8. A plot of the total energy E(total) as a function of bimeron separation a, for a layer 
separation of 0.7. This curve was obtained by extrapolating the curve fitted to the calculated 
values going upto a = 3.2 . 
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Table 1: The optimal bimeron separation (a), the bimeron local energy( Eiocai ) and meron 
pair energy {Eprev) from our previous work as a function of the layer separation d. The 

2 

unit of energy is ^ and the unit of length is / 



d 


a 


Elocal 


E 


0.5 


3.30 


.270 


.217 


0.6 


3.16 


.248 


.226 


0.7 


2.72 


.223 


.224 


0.8 


2.39 


.201 


.214 
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